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Abstract 

We compute the two-loop renormalization functions, in the RI' scheme, of local bilinear quark 
operators '^l'T^p, where F corresponds to the Vector, Axial- Vector and Tensor Dirac operators, in 
the lattice formulation of QCD. We consider both the flavor nonsinglet and singlet operators. 

We use the clover action for fermions and the Wilson action for gluons. Our results are given 
as a polynomial in csw: in terms of both the renormalized and bare coupling constant, in the 
renormalized Feynman gauge. 



Finally, we present our results in the MS scheme, for easier comparison with calculations in the 
continuum. The corresponding results, for fermions in an arbitrary representation, together with 
some special features of superficially divergent integrals, are included in the Appendices. 

PACS numbers: 11.15.Ha, 12.38.Gc, ll.lO.Gh, 12.38.Bx 
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I. INTRODUCTION 



Numerical simulations of QCD, formulated on the lattice, make use of a variety of com- 
posite operators, made out of quark fields. In particular, matrix elements and correlation 
functions of such operators, which include local and extended bilinears, as well as four- 
fermion operators, are computed in order to study hadronic properties in this context. A 
proper renormalization of these operators is essential for the extraction of physical results 
from the dimensionless quantities measured in numerical simulations. 

The present paper is the second in a series of papers regarding the calculation of renormal- 
ization functions of fermion bilinear operators to two loops in Lattice perturbation theory. 
The calculation of the scalar and pseudoscalar cases was carried out in Ref. In this work 
we study the renormalization function Zr of fermion bilinears O = tpTip on the lattice, where 
r = 75 7/i, l50'fiu {o'fj.u = 1/2 [ifiylu])- We consider both flavor singlet and nonsinglet 
operators. We employ the standard Wilson action for gluons and clover-improved Wilson 
fermions. The number of quark flavors Nf, the number of colors Nc and the clover coefficient 
cgw are kept as free parameters. One necessary ingredient for the renormalization of fermion 
bilinears is the 2-loop quark field renormalization, Z^, calculated in [l]. The one- loop ex- 
pression for the renormalization function Zg of the coupling constant is also necessary for 
expressing the results in terms of both the bare and the renormalized coupling constant. 

Our two-loop calculations have been performed in the bare and in the renormalized 
Feynman gauge. For the latter, we need the 1-loop renormalization functions Za and Za of 
the gauge parameter and gluon field respectively, as well as the one-loop expressions for Zr 
with an arbitrary value of the gauge parameter. 

The main results presented in this work are the following 2-loop bare Green's functions 
(amputated, one-particle irreducible (IPI)): 

• 2-pt function of the vector operator ip'j^ip : Sy(gaj^) 

• 2-pt function of the axial- vector operator ip^y^j^ip : Tij^y{qa^^) 

• 2-pt function of the tensor operator 4'j5cr^u4' '■ ^ril^h) 

(Oj^ : lattice spacing, q : external momentum). 

In general, one can use bare Green's functions to construct Z^' , the renormalization 
function for operator O, computed within a regularization X [X = L: lattice regularization; 
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X = DR: dimensional regularization) and renormalized in a scheme Y. We employ two 
widely used schemes to compute the various two-loop renormalization functions: 

• The RI' scheme: Zy^'' , Z^y^ , Z^'^^ 



The MS scheme: Z 



L,MS 
V ' 



7L,MS 



7L,MS 



For convenience, the results for Zq ' are given in terms of both the bare coupling constant 
Qo and the renormalized one: Qri', gjjs- Finally, as one of several checks on our results, we 
construct the 2- loop renormalized Green's functions in RI': ^^^'(g, /i) {O = V,AV,T), as 
well as their counterparts in MS: T,Q^{q,fL). The values of all these functions, computed 
on the lattice, coincide with values computed in dimensional regularization (we derive the 

n 

latter from the results of Ref. [2|). 

The present work, along with is the first two- loop computation of the renormalization 
of fermion bilinears on the lattice. One-loop computations of the same quantities exist for 
quite some time now (see, e.g., j^, Q], [sl and references therein). There have been made 
several attempts to estimate Zq non-perturbatively; recent results can be found in Refs. 



theory 
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Tl. A series of results have also been obtained using stochastic perturbation 
14l |. A related computation, regarding the fermion mass renormalization 



with staggered fermions can be found in 



15|- 



The paper is organized as follows: Section [ITl provides a formulation of the problem, as 
well as all necessary definitions of renormalization schemes and of the quantities to compute. 
Section IlIII describes our computational methods and results. Finally, in Section IIVI we 
discuss some salient features of our calculation, and comment on future extensions to the 
present work. 

Recently, there has been some interest in gauge theories with fermions in different rep- 



resentations 16 |o 

0, Q, Q, 
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the gauge group. Such theories are being studied in various contexts 



22|, e.g., supersymmetry [23|, phase transitions 



2^, and the AdS/QCD 



correspondence. It is relatively straightforward to generalize our results to an arbitrary 
representation; this is presented in Appendix A. Some special features of 2-, 3- and 4-index 
superficially divergent integrals are described in Appendix B. Finally, a detailed presentation 
of our calculation results on a per diagram basis, is provided in Appendix C. 
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II. FORMULATION OF THE PROBLEM 



Lattice action 



In the present work we employ the Wilson formulation of the QCD action on the lattice, 



with the addition of the clover (SW) 25| term for fermions. In standard notation, it reads: 



f ^ 



/ a:, At 

where : F^^ = ^ {Q^^ - Q^^) (2) 
and : 











p 


— 'z 9 i.Q^i' Quii) 








Ux,x+^Ux+^,x+^+uUx+^- 


-V, X- 


-V Ux-\-I/^ X 




~\~ Ux, x+u Ux-\-u, x+u—^ Ux-\-u- 


-fl, X' 


-fj, Ux — fM, X 




~^ Ux,x—^Ux—^,x—^—uUx—^- 


-V, x~ 


-V V X — V, X 




~l~ Ux^ X—V Ux—U, X~U+fi Ux—l/- 


-/I, x-\ 


-fi f^x+zx, X 



(3) 

Sg is the standard pure gluon action, made out of 1 x 1 plaquettes. The clover coefficient 
csw is treated here as a free parameter; r is the Wilson parameter (set to r = 1 henceforth); 
/ is a flavor index; a^j^u = (z/2)[7^, •y^]. Powers of the lattice spacing have been omitted 
and may be directly reinserted by dimensional counting. 

The "Lagrangian mass" rrio is a free parameter here. However, since we will be using 
mass independent renormalization schemes, all renormalization functions which we will be 
calculating, must be evaluated at vanishing renormalized mass, that is, when rrio is set equal 
to the critical value rricr'- rrio ^cr = rrii g1 + 0{gf). 

B. Definition of renormalized quantities 

As a prerequisite to our programme, we will use the renormalization functions, Za-, Zc, 
Z^, Zg and Z^, for the gluon, ghost and fermion fields (A^, c", ■?/'), and for the coupling 



constant g and gauge parameter a, respectively (for definitions of these quantities, see Ref. 
[l| ) ; we will also need the fermion mass counterterm . These quantities are all needed to 
one loop, except for which is required to two loops. The value of each Zq depends both 
on the regularization X and on the renormalization scheme Y employed, and thus should 
properly be denoted as Z'^'^ . 

Our one-loop results for the Vector and Axial- Vector operators, even though performed 
in a generic gauge, turn out to be independent of the gauge parameter. These results along 
with the one-loop expression for the Tensor operator, are in agreement with results found 
in the literature (see, e.g., Ref. sl). 



As mentioned before, we employ the RI' renormalization scheme 



26 
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which is 



more immediate for a lattice regularized theory. It is defined by imposing a set of nor- 



malization conditions on matrix elements at a scale yU, where (just as in the MS scheme) 



29|: 



/i = /i(47r/e^^)^/2 

where 7e is the Euler constant and ^ is the scale entering the bare coupling constant Qo = 
/i*^ Zg g when regularizing in Z) = 4 — 2e dimensions. 



C. Conversion to the MS scheme 

For easier comparison with calculations coming from the continuum, we need to express 
our results in the MS scheme. Each renormalization function on the lattice, Z^^^ , may be 
expressed as a power series in the renormalized coupling constant Qrii. For the purposes of 



our work the conversion of Qrii to MS is trivial since: 

9Bj' = ms + o{{msf) (5) 



The conversion of the gauge parameter apui to the MS scheme is given by [3 

^Ri' = Jlrv ^ms = ^MS I '^A^ms-. ^JTs) (6) 
where the conversion factor Ca may be calculated more easily in dimensional regularization 



(DR) 
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since the ratio of Z's appearing in Eq. ([H]) is necessarily regularization independent. 
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To one loop, the conversion factor Ca equals: 

^DR,RI' 2 

CaIq, a) = -^—= = 1 + / r(9Q;2 + 18a + 97) - 40Nf] (7) 

^ ^DR,MS 36(167r2) LV ' ' y c /j \ ) 



(Here, and throughout the rest of this work, both g and a are in the MS scheme, unless 
specified otherwise.) 

Once we have computed the renormalized Green's functions in the RI' scheme, we can 



construct their MS counterparts using the quark field conversion factor which, up to the 
required perturbative order, is given by: 



yL,MS r^DR,MS 



{8a^ + 5) CF + 14 Nf 
- (9q;^ - 24C(3) a + 52a- 24C(3) + 82) 7vJ(8) 

where cp — {N^ — l)/(2 Nc) is the quadratic Casimir operator in the fundamental represen- 
tation of the color group; ({x) is Riemann's zeta function. 



D. Renormalization of fermion bilinears 

The lattice operators Or = ipT ip must, in general, be renormalized in order to have finite 
matrix elements. We define renormalized operators by 

Of = Z^''''\a,i,)Oro (9) 

The fiavor singlet axial-vector Green's function receives additional contributions as com- 
pared to the nonsinglet case, while for the rest of the operators under study, singlet and 
nonsinglet Green's functions coincide. For the vector (V), axial- vector (AV) and tensor (T) 
operators, the renormalization functions Z^'^^ can be extracted through the corresponding 
bare 2-point functions 2^(50^) (amputated, IPI) on the lattice. Let us first express these 
bare Green's functions in the following way: 
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Siy(gaJ = ^5^^si^i.(gaJ +75^4y(9«J (10) 
S^(gaJ = 75 cr,^^? (ga J + i J^^''^^' ^^^^^ sg^ (ga J 

It is worth noting here that terms which break Lorentz invariance (but are compatible with 
hypercubic invariance), such as '^^{q^Y/q^-, turn out to be absent from all bare Green's 
functions; thus, the latter have the same Lorentz structure as in the continuum. Let us 
also point out that the presence of the 75 matrix in the tensor operator definition does 
not affect the bare Green's function on the lattice, in the RI' scheme. We have performed 
the calculation both with and without the inclusion of the 75 matrix, and we ended up 
with identical 2-point functions. Thus, for the purpose of converting our results to the MS 
scheme, we employed the conversion factors given in Ref. j^, where the definition of the 
tensor operator does not contain the 75 matrix. Furthermore, we expect that Yyj, {qa^) must 
vanish, since this is the case for the corresponding quantity coming from the continuum. 
Indeed, after performing the calculation on the lattice, it turns out that all contributions of 
this type vanish. 

Once all necessary Feynman diagrams contributing to the bare Green's functions pre- 
sented above are evaluated, one can obtain the renormalization functions for the three 
operators through the following conditions: 



where: 



lim 



lim 



lim 



7L,Rr ryL,RI' sn(l),L 



7L,RI' ryL,RI' sr(l),L 



7L,RI' ryL,RI' 



{qa^ 
{qa^ 
{qa^ 



7/. 

75 7m 

75 (^nu 



(11) 

(12) 
(13) 



^Avi(l(^L 

Ti^\qa^ 



1 + 0{gl 



sg^(ga,- 







(14) 



The conversion of the quantities Z-^ ' to the MS scheme is a straightforward procedure. 
In the case of the vector and tensor operators, the renormalization functions, Zy^^^ and 
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Z^'^^ , can be obtained by: 



Z^'^^ = Z^'^''/Crig,a) (15) 



where Cr{g,<y) are regularization independent conversion factors (F = V,T). These conver- 
sion factors have been calculated in dimensional regularization j^: 

r^L,Rr r^DR,RI' 

Cv{9. «) ^ ^ = = 1 + (16) 

Ct(5',") 

{2\^o? + 4320C(3) - 4815) cp - 626 Nf 



ryL,MS rvDRMS 

Zjrp Zjrp 

ryLjIS ~ ryDRMS 

Zjrp Zjrp 

, g" / 

1 H ^ Ci? a H ; r—r Cp 

167r2 ^ 216(167r2)2 



+ (I62a2 + 756a - 3024C(3) + 5987) iV, (17) 
Unlike the tensor operator, where the presence of the 75 matrix is irrelevant, the axial- 



vector bilinear {p\y = V'oTsT/x^o) requires special attention also in the MS scheme, due 
to the non-unique generalization of 75 to D dimensions. A practical definition of 75 for 
multiloop calculations, which is most commonly employed in dimensional regularization 



and does not suffer from inconsistencies is [Mfl 



75 = ^^e''^'^^"^'^*7.,7..7.3 7.4 , = 0, 1, 2, 3 (18) 

Of course, 75 as defined in Eq. f|T8l) does not anticommute (in D dimensions) with 7^, for 
/i > 4; an ultimate consequence of this fact is that Ward identities involving the axial- vector 
and pseudoscalar operators, renormalized in this way, are violated. 

To obtain a correctly normalized axial- vector operator [32] , O^y , one must introduce an 
extra finite factor, Z5, in addition to the usual renormalization function Z^y'^^^ (the latter 
only contains poles in e). We set: 

oW' = ZM oT = ZM z^^^^ o\y (19) 



For the definition of Z^ we must express the MS renormalized Green's functions Gy^ , 
G^y as well as the renormalized Green's function G^y (corresponding to O^y ), in a form 
similar to Eq. flTT]) : 

Cn^aJ = 7.G«^(gaJ + ^Gi,^)^(gaJ 



Gj'#(,aJ = 757, G2"'(?aJ + 7,^gS"'(?<'J (20) 



Z5 is then defined by the requirement that the renormahzed Green's functions Gy"^ (qa^) 

'AV 



and G^y^'^ coincide: 



G 



(1) Af5 



Z5 = ^ _ (21) 
'-'Ay 

Eq. f|T9l) is vahd for both the singlet and nonsinglet currents, provided of course, the appropri- 
ate choice for is used. Thus, we have two different expressions, Z| and corresponding 
to the singlet and nonsinglet axial-vector operator, respectively. They are gauge independent 
and differ only in the cpNf term; this is expected considering the fact that the additional 
Feynman diagrams contributing to the singlet axial operator have the insertion within the 
closed fermion loop. Both Z^ and Z^^ were evaluated in DR js^ and up to two loops they 
read: 

^7' „ ^ , / f.. 2 107 _ , 31 



^5%) - l-Ti^(4c.)+(^ [224--c.N.+ -c.N,j (22) 

Zr{9) = l-lf^(4c-)+(l|^ (22cl-'-^c,N^+lc,N,"j (23) 
^Av^^ can now be obtained by: 

ZAf^ = Z^Af'/ (Cv Z,) (24) 

where Z5 stands for Z| or Zg'^ (Eqs. fl22ti23p ). for the singlet or nonsinglet cases, respectively. 

Similarly, one can convert the RI' renormahzed Green's functions, Gy^\ to their MS 
counterparts, through: 

fiRI' r^RI' ^RI' 

^v ri ^ AV n n 7 t n n 

(j-y ^AV 

(In Eqs. ( l2Hl25i) it is understood that powers of qrp^ curi', implicit in RI' quantities, must 
also be converted to gjfs^ ^JIs^ respectively, using Eqs.(l5H6])). Note that the combination 
Gy Z^ appearing above yields the value of Gav = Z^v'^^ /Z^v'^^^ = Gy Z^. 
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III. COMPUTATION AND RESULTS 



The Feynman diagrams contributing to the bare Green's functions for the vector, axial- 
vector and tensor operators, Sy^yj,(g, a^), at 1- and 2-loop level, are shown in Figs. [Hand 
[21 respectively. For flavor singlet bilinears, there are 4 extra diagrams, shown in Fig. [31 
which contain the operator insertion inside a closed fermion loop. These diagrams give a 
nonzero contribution only in the axial-vector case. 

The evaluation and algebraic manipulation of Feynman diagrams, leading to a code for 
numerical loop integration, is performed automatically using our software for Lattice Per- 
turbation Theory, written in Mathematica. 

The most laborious aspect of the procedure is the extraction of the dependence on the 
external momentum q. This is a delicate task at two loops; for this purpose, we cast algebraic 
expressions (typically involving thousands of summands) into terms which can be naively 
Taylor expanded in q to the required order, plus a smaller set of terms containing superficial 
divergences and/or sub divergences. The latter can be evaluated by an extension of the 



method of Ref. 



33] to 2 loops; this entails analytical continuation to D > A dimensions. 



and splitting each expression into a UV-finite part (which can thus be calculated in the 



continuum, using the methods of Ref. 



34|). and a part which is polynomial in q. A primitive 



set of divergent lattice integral 
manner, can be found in Ref. 



s involving gluon propagators, which can be obtained in this 



35I . Due to the presence of at least one free Lorentz index in 
the definition of the operators (for the case of the tensor bilinear there are two such indices), 
it is possible to end up dealing with superficially divergent integrals with two, three or even 
four free Lorentz indices. In Appendix B, we provide a brief description of the manipulations 
performed to resolve such terms, based on the method described above. 

Some of the diagrams contributing to j^yrp[qa^) are infrared divergent when considered 
separately, and thus must be grouped together in order to give finite results. Such groups 



FIG. 1: One- loop diagram contributing to Zy, Zav and Zt- A wavy (solid) line represents gluons 
(fermions). A cross denotes the Dirac matrices 7^ (vector), 757^ (axial vector) and 750-^1/ (tensor). 



10 




19 20 

FIG. 2: Two-loop diagrams contributing to Zy, Zav and Zj-. Wavy (solid, dotted) lines represent 
gluons (fermions, ghosts). A solid box denotes a vertex from the measure part of the action; a 
solid circle is a mass counterterm; crosses denote the matrices jfj, (vector) , 757^^ (axial- vector) and 
750-/ii/ (tensor). 




12 3 4 



FIG. 3: Extra two-loop diagrams contributing to Zj^y, singlet- A cross denotes an insertion of a 
flavor singlet operator. Wavy (solid) lines represent gluons (fermions). 
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are formed by diagrams (3-7), (8-9), (10-11,19) in Fig. [2] and diagrams (1-2), (3-4) in Fig. 

m 

In Figs. [T] to [31 "mirror" diagrams (those in which the direction of the external fermion 
hne is reversed) should also be included. In most cases, these coincide trivially with the 
original diagrams; even in the remaining cases, they can be seen to give equal contribution, 
by invariance under charge conjugation. 

As mentioned before, all calculations should be performed at vanishing renormalized 
mass; this can be achieved by working with massless fermion propagators, provided an 
appropriate fermion mass counterterm is introduced (diagram 11 in Fig. [2]). 

All two-loop diagrams have been calculated in the bare Feynman gauge (do = 1)- One- 
loop diagrams have been calculated for generic values of Oo] this allows us to convert our 
two- loop results to the renormalized Feynman gauge [apii' = 1 or ajfg = 1). After perform- 
ing the calculation for the cases of the vector and axial-vector operator, we see that one-loop 
expressions for the renormalization functions do not depend on the gauge parameter. Espe- 
cially for the case of the vector operator, having in mind Eqs. (JME]) and Eq. (|T6l) . this fact 



causes the lattice results in the RI' and in the MS scheme to coincide. 

Numerical loop integration was carried out by our "integrator" program, a metacode 
written in Mathematica, for converting lengthy integrands into efficient Fortran code. Two- 
loop numerical integrals were evaluated as sums over finite lattices, of size up to L = 40; the 
results were then extrapolated to L ^ oo. Extrapolation is the only source of systematic 



error; this error can be estimated quite accurately (see, e.g. Ref. [36|), given that L- 
dependence of results can only span a restricted set of functional forms. 

A. One-loop results 

1-loop results for Z^'^^ are presented below in a generic gauge. As it turns out, only 
the tensor renormalization function depends on the gauge parameter, while for all other 
operators, one-loop expressions that emerge are gauge independent. The errors appearing 
in the expression for Z^'^^ , result from the L ^ oo extrapolation. 



«2 

ryL,Rr ^ I i/o 



- ln(a^/i^) + fto - 17.018079209(7) 
+3.91333261(4) cgw + 1.972295300(5) 4^ (26) 
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The corresponding expressions for Zy^^ , Z^y^ can be read off from Eqs. (130113 ip be- 
low. In all cases, one-loop results in the MS scheme {Zy'^^^ , Z^^^ , Z^'^^) present no 
dependence on the gauge parameter. 



B. Two-loop results 



In order to derive the expressions for the bare Green's functions S^^y^! ^'^^ must evaluate 
all Feynman diagrams presented in Figs. [T]l3l The extraction of Zy"^ , Z^y and Z^'^ is then 
straightforward via Eqs. (JIUIISD (for Y = RI'), and via Eqs. (UMIIl El (for Y = JdS). 
To this end, we need the following one-loop expression for Z^'"^ (note that = 1 to this 
order) : 



7L,RI' 



167r2 



+A^/ (-2.168501047(1) + 0.7969452308(4) cgw - 4.7126914428(1) 4^) 



+39.47841760436(1) cp + 1.94017130069(1) + O (g^) 



(27) 



To express our results in terms of the renormalized coupling constant, we also need the 
one- loop expression for 



7L,Rr 



_9l_ 
167r2 



+A^/ (0.5286949677(5) - 0.3984726154(2) cgw + 2.35634572140(7) 4^) 
-19.73920880218(1) cp - 3.54958342046(1) + O [gl) (28) 



Eqs. fi m [25]) are in agreement with older references (see, e.g., Ref. [37j). 

A final necessary ingredient is the two- loop expression for Z^'^^ , as required by Eqs. (iTTl - 
T3l) : this was calculated in Ref. l|], in the renormalized Feynman gauge aju' = 1, and is 
included here for completeness: 
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yL,Rr 



1 + 



+ 



167r2 

at 



Cf 



ln{alp,^) + 11.852404288(5) - 2.248868528(3) csw - 1.397267102(5) 4w 



■Cf 



1 2/ 2-2\ 

111 (a^/i ^ 
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(167r2)2 

+ ln(a^/i2) - 6.36317446(8) Nf + 0.79694523(2) Nf cgw 

-4.712691443(4) A^/4w 
+49.83082185(5) cf - 2.24886861(7) cf csw 
-1.39726705(1) cf 4^ + 29.03029398(4) N^^ 
+Nf - 7.838(2) + 1.153(1) cgw + 3.202(3) 4^ 

+6.2477(6) ciw + 4.0232(6) 4w) 
+CF (505.39(1) - 58.210(9) csw + 20.405(5) 4w 
+18.8431(8) 4^ + 4.2793(2) 4^) 
- 20.59(1) - 3.190(5) cgw - 23.107(6) 4w 
-5.7234(5) ciw - 0.7938(1) 4^) 

We present below Zy^^ , ^^^^ and Z^'^^ to two loops in the renormalized Feynman 



(29) 



gauge aRp = 1; we also present the MS analogues Z^^^ and Z^'^"'^ in the gauge — 
1 (as already mentioned, Zy^^ — Zy^^ ) . The bare Green's functions are relegated to 
Appendix C, where a per c/za^ram breakdown of the results is provided. It is a straightforward 
exercise to recover the total bare Green's functions from the corresponding Z's and the 
renormalized Green's functions. 

-.2 



LMS ■ 



1 + 



+ 



9o 



Cf 



9t 



(167r2)' 



Cf 



- 20.617798655(6) + 4.745564682(3) cgw + 0.543168028(5) c|w] 

Nf (^25.610(3) - 11.058(1) csw + 33.937(3) 4^ 

-13.5286(6) 4w - 1.2914(6) c^^) 
+CF ( - 539.78(1) - 223.57(2) csw - 104.116(5) 4^ 

-32.2623(8) c^w + 4.5575(3) c^w) 
+iVc - 51.59(1) + 18.543(5) cgw + 20.960(6) 4^ 

+2.5121(5) ciw + 0.1765(1) 4w)] (30) 
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yL,Rr 



= 1 + 



+ 



9l 



■Cf 



15.796283066(5) - 0.247827627(3) cgw + 2.251366176(5) 4 

Nf (^18.497(3) - 1.285(1) csw + 19.071(3) 4w 

+1.0333(6) 4w - 6.7549(6) c|w) 
+cf ( - 184.01(1) - 389.86(1) cgw - 166.738(6) 4^ 

+7.894(1) + 4.3201(3) 4w) 
+A^c ( - 21.62(1) - 33.652(5) csw + 26.636(6) 4mv 

+10.2186(5) 4^ + 1.4893(1) 4^) 



•sw 



(31) 



yL,MS 

■'AV 



= 1 + 



+ 



9l 
16^'^ 

(167r2)2 



11.796283066(5) - 0.247827627(3) csw + 2.251366176(5) 4^ 

+Ar^ (^14.045(3) + 1.903(1) cgw + 0.220(3) 4^ 

+1.0333(6) 4w - 6.7549(6) c^w) 
+CF - 95.28(1) - 390.85(1) cgw - 157.733(6) 4w 

+7.894(1) + 4.3201(3) 4w) 
+A^e (l8.67(l) - 33.652(5) cgw + 26.636(6) 4^ 

+10.2186(5) 4w + 1.4893(1) 4w) 



(32) 



' + 16^^^ 



+ 



9t 



(167r2)' 



- In(a^/i2) - 16.018079209(7) 
+3.91333261(4) csw + 1.972295300(5) 4w 
ln\allP) QcF - \Nf - Hat,) 

ln(aj//^) (^3.1685002(6) Nf - 0.79694524(6) A^/Cgw 

+4.71269143(3) iVjclw 
-13.96033835(7) cf - 3.9133325(1) cf cgw 
-1.97229535(2) cf4w - 25.04361149(6) A^^) 
+Nf (^16.923(6) - 8.399(2) csw + 18.711(3) 



+ 



-sw 
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-10.8351(8) 4w - 5.1253(6) 
+CF ( - 868.0(1) + 551.6(2) cgw + 63.9(1) 4^ 

-79.49(1) 4^- 12.586(1) c^w) 
+N, - 15.76(8) + 27.6(1) csw + 38.2(1) 

+7.021(8) c|w + 1-6653(9) 4w 



(33) 



^2 

167r^ 



- ln(a^/i^) - 17.018079209(7) 
+3.91333261(4) cgw + 1.972295300(5) c^^ 



ln^«/.^)(-iiV^ + icp + ^iV, 



+ ln(a^/i2) (^2.5018336(6) Nf - 0.79694524(6) Nj cgw 

+4.71269143(3) A^/4w 
-12.96033835(7) cp - 3.9133325(1) cp csw 
-1.97229535(2) c^c^w " 21.37694482(6) A/", 
+Nf (^21.989(6) - 9.196(2) cgw + 23.424(3) c^^v 

-10.8351(8) 4w - 5.1253(6) 4w) 
+CF ( - 893.2(1) + 547.7(2) cgw + 61.9(1) 

-79.49(1) 4^- 12.586(1) 4w) 
+N, - 41.44(8) + 27.6(1) cgw + 38.2(1) ci 



+7.021(8) 4w + 1-6653(9) c 



sw 



(34) 



All expressions reported thus far for Zy, Zav and Zt refer to flavor nonsinglet operators. 
In the case of Zy and Zt, all diagrams of Fig. [3] vanish, so that singlet and nonsinglet 
results coincide, just as in dimensional regularization. For Zj^y on the other hand, the 
above diagrams give an additional contribution: 

4 

Igiet = ^Ay + ^^^c^AT^ (-6 Ma^/x^)- 2.0491(5) + 15.0315(6) Csw 

+5.0090(2) 4w - 2.11016(5) 4^ - 0.04329(2) 4w) (35) 



The same extra contribution applies also to the MS scheme. 

For the sake of completeness, and as an additional check on our results, we compute the 
renormalized Green's functions (for vanishing renormalized mass). Since the bare Green's 
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functions have two contributions of different structure (as defined in Eq. ( JTOj) . see also Eq. 
fl20|) ). we derive the renormahzed expressions for these contributions separately: 



G 
G 



G 



(i) RV 
V 

(i) RI' 
AV 

(i) RI' 



yL,Rr yL,Rr ^{i),L 

yL,RI' yL,Rr sp(.i),L 
^AV ^AV 



(36) 
(37) 
(38) 



where i = 1,2. Similarly for MS, taking into account Eq. ( ]25ll . 

Since these functions are regularization independent, they can be calculated also using, 
e.g., dimensional regularization. We have computed G^ , G^y and G'^^ in both ways: either 
starting from our Eqs. (12611341) or using renormalization functions from dimensional regular- 
ization [2I. In all cases the two ways are in complete agreement. We obtain: 

-.2 



G?^^^'(?) = l + ^c^ln(/.V?^) 
167r^ 



+ 



9%! 



G 



(2) RI' 
V 



(?) 



+ 



(167r 



9ri 
167r 

9ri 



2^2 



cf 



(167r 



ln^(^V?') 

+ ln(/iV?') 
-CF (-21n(/iVg')) 
Hf/q') (• 
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Nf cf + 

9^2 1 



251 



2^2 



Cf 



-2cF-^N,) + ^Nf + 3cF 



251 



AT, 



(39) 



(40) 



The vector renormahzed Green's function in the RI' scheme coincides with the corre- 
sponding axial- vector expression, and thus Eqs. fl39tl40|) also hold for the axial- vector case: 
^Aviq) — C!v\q)^ ^Aviq) — C^y\q). Of course, even though the MS expression for the 
vector renormalization function, Zy'^^ , coincides with the RI' expression, that is not the 



case for the renormahzed MS Green's function, due to G^ appearing in Eq. ( !25|) . This 
factor results in the following quantities: 



G 



(1) MS 
V 



(?) 



167r2 

'MS 



CF (ln(/.V?') + 1) 



q— 



'167r 



2^2 



Cf 



\^\^Vq')(\cF + K 



+ ln(^V?') 
7 
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Nf 



1 251 
-cf^ Nc 

2 18 



+ 



/ 



■Cf 



143 



6C(3) 



(41) 
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2^2 



Cf 



(42) 



Furthermore, the axial-vector renormahzed 2-point functions in the MS scheme differ 
from Eqs. (14 11142 p . due to the finite conversion factor Z^^; they read: 

-.2 



9- 



167r2 

•,4 



+ 



9- 



MS 



(167r 



2^2 



Cf 



w 2 / 2N 19 7 251 \ 

+ ln(/iVg2) (^-_ iV^ + - + _ AT, j 



(43) 



+ 



^2CF 

Cf 



167r2 



y y 



(44) 



(167r2)2 

If one considers the singlet axial- vector current, then there exists an extra contribution to 
the expressions above: 

^(l),MS / ^ _ ^(1)M5/ ^ 9^ 

singlet W J - '-^AV + /iR-,r2^2 



G 



{2),MS 
AV, singlet 



WS, ^ . 9^ 



-6 In(^Vg^) cfNj - - cfNj 



{<!) = gT\<^ + 



(167r 



2A2 



/J 



(45) 
(46) 



For the RI' scheme, similar relations hold, the only difference being the absence of the 



factors Z|, Zg^; we obtain: 



y(l) RI' 



9%P 



G 



{2),RI' 



9%r 



-4 CFNf) 



'AV, singlet V'iy ^ AV y-i/ • (^lQ^2y 

Finally, for the tensor renormalized Green's function, we obtain: 



(47) 
(48) 



G'J)^^'(g) = l + 



(167r2)2 



Cf 



7 



+ \nm'/q'){--Nf + 8cF--N, 



(49) 
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Just as was expected from dimensional regularization, G^^^^ (q) = 0. The corresponding 



quantity in the MS scheme reads: 

(1677 



.iiV,.(f-20C(3))c..(-^.SC(3)) 



(50) 



In Figs, a dSlD, dZHH]) we plot Z^'"^'', {Z^f\ Z^f^) and (Z^'^^', Z^'^^), respectively, 
as a function of csw- Values of the clover parameter used in simulations lie within the typical 
range < csw ^ 2. For definiteness, we have set = 3, ft = and Po = '^Nc/qI = 6.0. 
Our results up to two loops for each Z are shown for both Nf = and Nf = 2, and compared 
to the corresponding one-loop results. Furthermore, in the axial-vector case, we also present 
the two-loop result for the flavor singlet operator, for Nf = 2. 

In Fig. [9] we present, on the same plot, the values of Zy^^ , Z^y^ , ^^if^mc/iei ^^"^ Z^'^^ 
up to 2 loops, versus csw- We have chosen Nc = 3, = Nf = 2 and Po = 5.3. The 



corresponding results in the MS scheme are plotted in Fig. [101 
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IV. DISCUSSION 



In this paper we have reported results regarding the Vector, Axial- Vector and Tensor 
fermion bilinear operators. This work, along with a previously published paper |l| regarding 
the Scalar and Pseudoscalar operators, provide a complete two loop calculation for the renor- 
malization functions for local fermion bilinears, considering both the singlet and nonsinglet 
cases. The two- loop wave function renormalization constant, Z^, which is a prerequisite for 
our calculation, was presented in Ref. [l| (the reader should also refer to this paper for any 
necessary notation not included in the present sequel paper). 

It is clear from Figs. H] to [S] that the two- loop renormalization functions differ significantly 
from 1-loop values; this difference must then be properly taken into account in reducing 
systematic error in MC simulations. At the same time, 2-loop contributions are typically 
smaller than 1-loop contributions, especially for csw ^ 1; indicating that the (asymptotic) 
perturbative series are under control. 

The results are presented as a function of the clover parameter, where the values of csw 
lie within the standard range < Csw < 2. Optimal values for Csw, which have been 
estimated both non-perturbatively [38| and perturbatively (to 1-loop) |25|, lie within this 
range. A breakdown of our results on a per diagram basis is presented in Appendix C, for 
completeness. 

As already mentioned, we take into account both singlet and nonsinglet operators. After 
evaluating all Feynman diagrams involved, we found that, for the Vector and Tensor oper- 
ators, singlet renormalization functions coincide with nonsinglet ones. On the other hand, 
the Axial- Vector operator receives an additional contribution in the flavor singlet case. 

The numerical integrations over loop momenta were executed on a Pentium IV cluster; 
they required the equivalent of 60 months on a single CPU. 

A possible extension to the present calculation is the renormalization of more extended 
operators, with the same continuum limit as we have considered here. A standard basis of 
higher dimension operators, with the same quantum numbers as the local bilinears which 
we have considered, can be found e.g. in Ref. Such operators are frequently used 

to reduce 0(0^) effects. A number of additional Feynman diagrams must be introduced, 
since the vertices coming from these operators may also contain gluon lines. However, the 
additional integrals resulting after the contractions will be free of superficial divergences. 
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leading to a less cumbersome computation, despite an increase in the size of the integrals. 
Further directions regard higher dimensional operators, such as ^ D'^ ■ ■ ■ D'^ r\E', which enter 
structure function calculations, and 4-fermion operators. 

Finally, our computation can be easily extended to improved lattice actions. With re- 
gard to improved fermion actions, such as those containing twisted mass terms 39|] or 
Osterwalder-Seiler terms {4^, our results remain unchanged, since they pertain to mass- 
independent schemes. Improving the gluon action, on the other hand, is more CPU con- 
suming, but conceptually straightforward: Splitting (in iterative fashion) the Symanzik 
propagator into a Wilson gluon propagator plus the remainder, leads to the same bare 
Green's functions as the ones presented in this paper, with the addition of superficially 
convergent terms, which can be more easily manipulated. Based on our experience with 
other similar calculations, the algebraic expressions for the integrands will grow roughly by 
a factor of 5; furthermore, the gluon propagator must now be inverted numerically for each 
value of the momentum, leading to an additional factor of <2 in CPU time. Finally, if one 
wishes to employ more than one set of values for the Symanzik coefficients, CPU time for 
numerical integration will increase almost proportionately. 
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APPENDIX A: FERMIONS IN AN ARBITRARY REPRESENTATION 



The results presented up to this point regarded renormahzation constants of various 
fermion bihnear operators constructed with fermions in the fundamental representation of 
the gauge group. Our results were expressed in terms of the clover parameter, csw; the 
number of fermions, Nf, the number of colors, Nc, and the quadratic Casimir operator in 
the fundamental representation, cp- 

Recently there has been interest in theories with fermions in other representations; some 



preliminary non-perturbative calculations have also appeared (see e.g. 4l|, 42|). In this 
Appendix we describe the conventions we use in our work, regarding the generators of the 
algebra, and we then express our findings in an arbitrary representation. 

Our results for Zy, Zav, Zt, Eqs. fl30 l [311 [331) . can be easily generalized to an action with 
Wilson/clover fermions in an arbitrary representation R, of dimensionality rf/j. 

In this case, the gluon part of the action remains the same, while all link variables 
appearing in the fermion part of the action assume the form: 

t/x,x+M = exp(z(7o^;:(a;)n U,,,+^ = expitgoA';^ix)T^) (Al) 

Using standard notation and conventions, the generators in the fundamental representa- 
tion satisfy: 

AT2 _ 1 1 

[T^ T^] = i r''" T", J2 ^"^^ = 1 ■ Cf = 1 ■ — , tr(T°r'') = 5^'' tp = S^'' - (A2) 

a 

In the representation R we have: 

a 

where: tR= {(IrCr) / {N^ - I). 

For the 1-loop quantities, Eqs. (127 1 [28l) . converting to the representation i? is a straight- 
forward substitution: 

Nf-^Nf {tn/tp) = Nr{2 tn) (A4) 
and, in addition, for Eq. f[2^ : 

Cf — ' Cr (A5) 

Aside from these changes, all algebraic expressions (and the numerical coefficients resulting 
from loop integrations) remain the same. 
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A similar reasoning applies to the 2-loop quantities in Eqs. (ISOj [STj [33]) : For most 
diagrams, once their value is expressed as a linear combination of cjp, cpNc and cpNf, it 
suffices to apply substitutions (lA4p and (lASp . The only exceptions are diagrams containing 
a gluon tadpole [diagram 3 of Fig. [21 1-loop diagrams, when expressed in terms of ar/', 
(^Rr by means of Zg, Za]'- In these cases, only one power of cp should be changed to cr; 
a possible additional power of originates from the gluon tadpole and should stay as 
is. This peculiarity implies that, in order to perform the substitutions as described above, 
one must start from the per diagram breakdown of 2-loop results. To avoid a lengthy 
presentation, we apply, instead, substitutions ( lA4p and ( 1A5I) indiscriminately on Eqs. (l30l 
EHinS]); consequently, we must then add a correction term, as follows: 

\r ~ \cF-*CR, Nf-*2NitR 

' cr {cr - cf) ■ [813.9580654(2) - 187.3473843(1) cgw 



(167r2)2 

-21.4434142(2) (A6) 



^AV \r ~ ^AV lcF-»Cfl, Nf^2Nftit 



+ CR {cr - Cf) ■ [623.6122595(2) + 9.7838425(1) cgw 



-88.8803741(2) c^w] (A7) 



ryL^RI' I fyL^RV I 

\r - \cF^cn, Nf-,2Nftii 



0^ 

Cr {cr - Cf) ■ [47r2 ln(af /i^) + 632.3684202(3) 



(167r2)2 

-154.492179(2) cgw - 77.8630975(2) (A8) 

Actually, the reader could arrive at these results without knowledge of the per diagram 
breakdown, by virtue of the following fact: All "exceptional" powers of Ci;' cancel out of 
Zy^^ , Z^y^ , Z^'^^ , if these are expressed in terms of the renormalized coupling constant 
aRi'. Thus, one may: 

• Express Eqs. ( !30l [311 \33\) in terms of qrj/ by means of = {Zg'^^') Qri', with Z^'^^' 
in the fundamental representation (Eq. ( [281) ) 

• Apply substitutions ( 1A4I) . ( ]A5[) throughout 

• If desired, reexpress everything in terms of Qo (using 

(^Z^^Ri'yi from Eq. (l28D, with 

Nf 2NftR and cf a.s is) 
No correction terms are necessary in this procedure. 
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APPENDIX B: MANIPULATION OF SUPERFICIALLY DIVERGENT AND 
SUBDIVERGENT TERMS 

In the case of the Scalar and Pseudoscalar bihnears ^ all superficially divergent terms 
involved at most one free Lorentz index, but when one considers other bilinear operators 
different structures may arise. For example, in the case of the vector and axial-vector oper- 
ators, two-index integrals may arise and, of course, when working with the tensor operator 
three- or even four-index integrals appear during the manipulation of superficially divergent 
terms. The most laborious aspect of such an evaluation is to extract the explicit depen- 
dence of the bare matrix element on the external momentum, by expressing the superficially 
divergent parts in terms of known primitive divergent integrals. 

We will focus on an arbitrary four-index integrand emerging, for example, from a 
"diamond" -like diagram with a F = 75 0"^,^ insertion. Taking into account the symmetries of 
this object, we can deduce all possible tensor structures which may appear, as a linear com- 
bination, in the result for the corresponding integral. Such structures are certainly tensors 
under the hypercubic group, but not necessarily so under the full 5*0(4) Euclidean rotation 
group: Terms such as S^^p'^ or q^/{q^Y flight be present^ (g: external momentum), and if 
so they might spoil the renormalizability and/or the Lorentz invariance of the theory. We 
must show that in all cases, such terms are absent. 

Let us begin by taking as an example an algebraic expression which contains both su- 
perficial and sub divergencies; this example serves as a prototype for all the cases we have 
encountered. Such an expression may arise from a "diamond" -like diagram with the insertion 





,4 ,4; 1/ op ,o(T 

a pa k k p p k 



(Bl) 



^ (p2)2 p + q k — p k + q 



where q is the external momentum and 






(B2) 



No summation is implied over the indices /x, z/, p, a. 



1, fi^iy = p^a; S^'^p'' 



0, otherwise; g"* = J2fiiQ'^)^- 
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From simple ultraviolet power counting on the term above, one can realize that the super- 
ficial degree of divergence is —8 and the degree of divergence in each of the two loops is —6 
and —4. Thus, this term is not only superficially divergent but also contains a subdivergence 
in the right loop. All divergences are resolved by using a BPHZ procedure. The potential IR 
divergences, which may arise in intermediate steps, necessitate working in D > 4 dimensions 
as in 331]. Performing a BPHZ subtraction for the right loop, we split the integral into two 



parts: 



-'sub WJ — 



(2^1") {jp'Y p + q k — p kP' k + q 

d^p p^ p^ f d^k /e^/c*^ 



(27r)^ (^2)2 J (27r)^ (P)3 J ^^^^ 

The last term in Eq. flB3p is a separable integral. The integral over momentum p is a 
standard primitively divergent integral (see, e.g., [ssl), whose value contains only Lorentz 
invariant structures. The integral over momentum k does not depend on the external mo- 
mentum g, and gives nonzero result only when the indices ^ and a are in the same direction. 
Thus, this term assumes the following structural form: 

In the remaining part of the original expression, we must still perform an extra subtraction, 
to cure the superficial divergence: 

/sub(g) = [/sub(g) - /sub(o)] + /sub(o) (B6) 

According to the BPHZ procedure, the quantity [/sub(Q') — -^sub(O)] is now UV-finite, and thus 
it equals the corresponding continuum expression. Consequently, once again, only Lorentz 
invariant structures arise: 

[/sub(g) - /sub(O)] — ^ a' "I^J^l + ^ {6^" q" + 6''-qPqf^ + 5>'P + qf^ q-) 

+c' 6^"" + d' 6"" + e' {6^"'6'''' + 6'"'6^'P) + f'6^'''6''P (B7) 
The last part of the integral I^^'"^{q) (last term in Eq. (]B6I) ) equals: 

J (2vr)(2D) (p2)3 (p)2 [^^^ p j ^^'^ 
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This g-independent integral could give rise to a structural form of the type 8^'^'^" ^ which 
would spoil Lorentz invariance; however, this problem is avoided since the indices /x, u 
in Jt^'^P^ actually originate from the insertion F = 75cr^jy. As a consequence, only the 
combination Jt^'^P'^ — ji^pp^ appears in the Feynman diagram, and no 3^'^''" contribution 
survives. Thus, we are led to: 



[^^"(0) - Cr'lO)] ^ a {5^P5'"' - d'Pd^"') (B9) 

Having completed the whole procedure for the integral shown in Eq. (IBip . we conclude that 
the only functional form that a four-index object (with the symmetries described above) can 
have reads: 



q2 



+ ^^{S^'^q'^qf -d^^q^qf^ (BIO) 



q^ 

We emphasize again that, even though the above expression would be obvious in a continuum 
regularization, it is not so on the lattice, where one could have ended up with terms breaking 
Lorentz invariance. 

Using similar considerations, we can prove that the two- and three-index expressions 
which appear in our calculation, will take the same structural form as in the continuum; 
i.e., they will be free of Lorentz non invariant contributions, which could be present a priori, 
such as 5^"' (g^)Vg^ 5^"'P, etc. 

Once we establish the structural form of the two-, three- and four-index integrals, we 
must compute the coefficients multiplying each tensor structure, such as the coefficients 
A{q'^), B{q^), C{q^) of Eq. fIBlOp . We illustrate the procedure by taking as an example the 
following two-index integral: 

/'"'(.) = / J^'fJlf' (BID 

J l^vrj (^2)2 p + g k + q 

Along the same lines of reasoning as above, we conclude that the lattice integral I'^'^ is 
of the same form as its continuum counterpart: 

I^'^iq) = ASf"" + (B12) 

q^ 
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The problem is now reduced to evaluating the coefficients A and B. Upon contracting the 
integral shown in Eq. flB12p with 5^^^ and q^qu, we get: 



/i = J^S'""!^" = DA + B 
h = J2q^q,I'"' = Aq'' + Bq'' 



(B13) 
(B14) 



where D is the number of dimensions (on the lattice, D = 4). Once we evaluate the integrals 
/i and I2, we are able to determine the quantities A and B through the following relations: 



(B15) 



Let us proceed with the evaluation of the integral Ji by contracting 6'^'^ with Eq. fIBlip . 
keeping only terms of order 0{q^). At first, we aim to reduce the number of propagators 
appearing in the denominator by employing the following property, which is valid on the 
lattice: 



k.p = - 



(BI6) 



Omitting some intermediate steps, the resulting expression can be written as follows: 



where 



h=Ia + h + Ic + h 



1 [d^pd^k {l^-p)Eoik'?iP'? 



(B17) 



{jPY p + q k—p kP' k + q 



k 



.op 



d^k 



d'^p 



op 

p 



1 fd^pd^k 



jO o 

k ■ p 



,2^2 



p + g 



k'^ k + q p^ p + q k — p 
d^pd'^k k - p If d'^pd'^k 



(27r)8 ^^2)2 kTt k^p 2 



{k-q) ■ jp-q) 

V J ip — q \ k'^ k — p 



(B18) 



(B19) 



(B20) 



(B21) 



As can be seen from the expressions above, we have managed to reduce diamond-like 
expressions into simpler integrals. Integral la is IR convergent: one can set g = and carry 
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out the integration numerically. Integral /;,, shown in Eq. (lB19p . is the product of two 1 



loop integrals: The first one is well known and tabulated in 



35|, whereas the second, being 



UV-finite, equals its continuum analogue and can be solved by using the following formula 



found in Ref. 



2\2-e-a-l3 



where £= (4-L))/2, Dp e 
/?, n, a) = (47r)' 



(47r2 



(T>0 



1 



Vn{p) 



(B22) 



J p=q 



; d"^ /dpf,dp^, and 
T{a + (3 -a -2 + e) 



B{2-e-a + n-a,2 



/3 



r(a) is the Gamma function and B{a, 13) = T{a) T{j3)/T{a+j3) is the Beta function. Vnip) is 
an arbitrary homogeneous polynomial inp: Vn{^p) = A"P„(p). All UV-convergent integrals 
in our calculation can be treated using Eq. ( ]B22p . for various values of a, P, Vnip) (and also 
Eq. (3.4) of [34] for diamond diagrams). Integral Jc of Eq. (lB20p can be treated with one 
further application of Eq. (]B16p . leading to either IR convergent contributions or tabulated 
integrals (such as the scalar "eye" integral, Eq. (C.5) of |35|). Regarding the integral 1^, 
employing integration by parts, we find: 



^ 1 fd^pd^k J2pCOs{k-q) 



(2vr) 



p^ p — q k'^ k — p 



qP 



{k-qy 



d^p d^k 
(27r) p2 p _ q 1^2 _ p 



(B23) 



Simple trigonometry on the first term in Eq. ( 1B23I) : Yip cos{k — q)^ = 4 — k'^/2 + 0{q) leads 
to the scalar eye diagram plus a simple, separable integral. The second term in Eq. ( ]B23I) is 
the derivative of a vector eye diagram, and can be resolved in two steps: (i) The integrand, 
[k —q)^/ {p^ p — q k'^ k — p ), after simple trigonometry and use of the symmetry k ^ p — k, 
may be expressed in terms of superficially convergent and/or known divergent integrals, plus 

2 ^ — ^ 2 

an integrand of the form p'^ / [jp' p — q k"^ k — p ); (ii) the latter, upon contraction with q^ 
and use of Eq. ( 1B16I) . is expressed completely in terms of known integrals. 

Some superficially convergent integrals, which contain sub divergences, often appear in 
various stages of our calculation. A simple prototype example is: 



p^ p 



■^2 - 2 

q k"^ k — p 



which is logarithmically divergent for g — > 0. In such subtraction of the form: 

1 1 / 1 1 



k — p 



1 

k'^ 



,k — p 
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leads to known separable integrals, plus terms in which one can set g = without appearance 
of divergences. 

In conclusion, using the steps which we outlined above, we have managed to evaluate 
the integral Ji of Eq. fIBlSp . by reducing diamond-like structures into simpler ones, leading 
to expressions containing UV-finite integrals and standard primitively divergent integrals 
whose values are known. Using similar considerations, one can also evaluate the integral I2 
of Eq. flB14l) . which is needed up to order 0{q^). With the evaluation of these two integrals, 
we can fully determine the coefficients A and B of Eq. (IBlSp . leading to the calculation of 
the original two-index integral. Let us point out that, throughout the whole procedure, the 
necessity to work m D ^ A dimensions does not emerge (It was only necessary in order to 
carry out demonstrations, such as Eqs. flB3tlB9p . leading to the conclusion: Eq. flBlOp ). 
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APPENDIX C: PER DIAGRAM RESULTS 



In this appendix we present our perturbative results for the bare Green's functions, 
Sr(5'o, (where T = V, AV, T), on a per diagram basis. Our results are expressed 

in terms of the bare coupling constant, Qo, the lattice spacing a^, the external momentum, 
q and the clover parameter csw- For the sake of simplicity we have set A^^c = 3 in two- loop 
expressions; at one-loop level the number of colors is left unspecified and the bare gauge 
parameter, cto, niay take arbitrary values. In all cases, the number of flavors, Nf, can take 
arbitrary values. 

Only one Feynman diagram, shown in Fig. [H contributes to one-loop expressions. Our 
corresponding results for the three operators read: 



^v,i-ioop{go,a^^q) — gl 



{Nl - 1) 



7m 



327r2 



log (aj q^) + 0.0151728775487(3) 



On 



+0.01258087658(1) - 0.007905256548(1) cgw 
-1-0.0027043227859(1) c^w' 



^Av,i-ioop{go-,ai^q) — go 



2 {Nl - 1) 



AT, 





%4 


(- 1 


+ 




V 167r2 








75 


7a. ( 


327r2 



(CI) 



log [al g2) + 0.0151728775487(3) 



Olo 



-0.002685425493(8) + 0.007905256548(1) cgw 
-0.0027043227859(1) 4w 



%4 



1 



^T, 1 



q^ \ 167r^ 
1 — a 



(C2) 



32n' 



- log {al q^) + 0.01200659055973(9) 



+0.00118313174(2) - 0.0052701710(1) csw 



-0.001820704300(1) c; 



sw 



(C3) 



The contribution to the bare Green's functions from the £-th two-loop diagram, can be 
written in the folowing form: 

V 1152^^ 



^y,L/oop(^o, g) = ^o^; 



7a. 



E4w( TT^ri log^a^) + v2 log + 



4" 



'SW 



2887r^ 



log (ajg2) + 



4,j 



(C4) 
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^AV,2-ioopi9o,a^q) = gtN}l5 



W {ay) +wi^} log {alq') 



3,i 



'SW 



w 



2887r^ 



log {ay) + 



w 



(C5) 



2-«oop(^0' «L = ^/ 75 (X, 



X 



'SW 



11527r4 



log^ {ay) + log {ay) + 



X 



3,i 



(C6) 



where the index i runs over all contributing two-loop diagrams. The dependence on csw 
is polynomial of degree up to 4 {i = 0, ■ ■ ■ , 4). The number of flavors, Nf , is raised to 
the power k where, of course, k = 1 only for diagrams 8 and 9 of Fig. [21 since they are 
the only diagrams containing a closed fermion loop; the remaining diagrams have k = 0. 
The coefficients w^^^ and x^^-* are numerical constants obtained upon evaluating each 
two- loop Feynman diagram. In Tables [fl UTl IIIUlTVl and IVTlVIl we present our results for v^^-*, 
w^^^ and x^^\ respectively, with accuracy up to 10 decimal places. 

In the case of the singlet Axial- Vector operator, it turns out that the only additional dia- 



gram contributing to S 



AV, singlet 



a^q) is diagram 4 of Fig. [3l It is a straight forward exercise 



to recover the bare matrix element; starting from Z^yl-,^^^^^ (Eq. fl35|) and Z^j^'^^ (Eq. fl29l) 
one can employ the RI' renormalization condition, Eq. 
matrix element. 



L,RI' 

to extract the corresponding 
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TABLE I: Contribution of two- loop diagrams to "^1^2-1000 = 1~" 



f 


/ 






(f) 


(f) 










0.0010901413 


-0.00955549(5) 


0.0021802826 


1 


1 








0.0089908(2) 







2 








-0.00522028(9) 
















0.00230566(3) 





2 


1 
















2 








-0.000041360(2) 













-0.0004010149 


0.01183022(2) 


-0.0008020299 


3-7 


1 








-0.00738844(3) 







2 








0.00260794(1) 













0.0000534687 


-0.0003947(1) 


0.0001069374 




1 








0.00032545(2) 





8-9 


2 








-0.00077830(1) 







3 








0.000389303(9) 







4 








-0.000146063(2) 













-0.0017442261 


0.0152888292 


-0.0034884521 


10-11 


1 








-0.0043548025 







2 








0.0014897419 






















12 


1 








0.000063005(3) 







2 








-0.000061816(2) 













0.0000390267 


-0.000161967(5) 


0.0000780534 


13 


1 





0.0001255891 


-0.00099763(1) 


0.0002511781 




2 





-0.0002000264 


0.00009461(1) 


-0.0004000527 




3 








-0.000046762(1) 
















0.00045044(3) 





14 


1 








-0.00147021(3) 







2 








-0.000043209(1) 







3 








0.000057661(2) 
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TABLE II: Contribution of two-loop diagrams to ^v2-ioop ('^ ~ 





/ 


vE 










27 


-0.0020904154 


0.0061629(1) 


-0.00418077(4) 


15 


1 





-0.0005507191 


-0.0010101(1) 


-0.0011014381 




2 





0.0001759227 


0.00026474(6) 


0.0003518455 




3 








0.000424491(6) 













0.0000390267 


-0.000342083(4) 


0.0000780534 


16 


1 





0.0001255891 


0.001157216(6) 


0.0002511781 




2 





-0.0002000264 


0.00178325(2) 


-0.0004000527 




3 








-0.000076152(1) 













-0.0000356458 


-0.000049925(1) 


0.00001874(2) 




1 








0.000198131(2) 





17 


2 








-0.000064731(4) 







3 








-0.0000161812(2) 







4 








0.0000032674 










-1 


0.0000558195 


0.000098634(7) 


0.00002159(5) 




1 





-0.0000211097 


-0.000057416(3) 


-0.0000422194 


18 


2 





0.0000249033 


-0.000179175(1) 


0.0000498067 




3 








0.0000335005(3) 







4 








0.0000077358(2) 










-8 


0.0002681483 


-0.0086175(2) 


-0.0000340358 




1 





0.0003206506 


0.0144022(6) 


0.0006413012 


19 


2 





0.0001992266 


0.00454038(9) 


0.0003984533 




3 








0.00128588(2) 







4 








-0.00047870(1) 










4 


-0.0008031274 


0.00453973(4) 


-0.00103589(2) 




1 





0.0001779933 


-0.00234305(6) 


0.0003559866 


20 


2 





-0.0000608900 


0.00090217(4) 


-0.0001217801 




3 








0.000065108(4) 







4 








0.0000428543(8) 
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TABLE III: Contribution of two-loop diagrams to ^AV2-ioop = 1 





?' 


(^) 
"'1.' 






"'1., 










0.0010901413 


-0.004299380(2) 


0.0021802826 


1 


1 








-0.0089908(2) 







2 








0.00522028(9) 
















-0.00010270(1) 





2 


1 
















2 








0.000041360(2) 













-0.0004010149 


-0.00177094(2) 


-0.0008020299 


3-7 


1 








0.00738844(3) 







2 








-0.00260794(1) 













0.0000534687 


-0.0000143(1) 


0.0001069374 




1 








-0.00019712(2) 





8-9 


2 








0.00001658(1) 







3 








-0.000389303(9) 







4 








0.000146063(2) 













-0.0017442261 


0.0068790161 


-0.0034884521 


10-11 


1 








0.0043548025 







2 








-0.0014897419 






















12 


1 








-0.000063005(3) 







2 








0.000061816(2) 













0.0000390267 


0.00017234(2) 


0.0000780534 


13 


1 





0.0001255891 


0.00172937(5) 


0.0002511781 




2 





-0.0002000264 


-0.00147678(5) 


-0.0004000527 




3 








0.000046762(1) 
















0.000203052(3) 





14 


1 








0.00057021(3) 







2 








0.00010067(1) 







3 








-0.000057664(2) 
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TABLE IV: Contribution of two-loop diagrams to '^AV2-ioop = 15 





i 












27 


-0.0020904154 


0.0019129(1) 


-0.00418077(4) 


15 


1 





-0.0005507191 


0.0048214(1) 


-0.0011014381 




2 





0.0001759227 


-0.00100389(6) 


0.0003518455 




3 








-0.000424491(6) 













0.0000390267 


-0.0001539166(3) 


0.0000780534 


16 


1 





0.0001255891 


0.00125957(2) 


0.0002511781 




2 





-0.0002000264 


0.001151401(6) 


-0.0004000527 




3 








0.000076152(1) 













-0.0000356458 


-0.000281627(7) 


0.0005068330 




1 








0.000060887(2) 





17 


2 








0.000050595(4) 







3 








-0.000006587(2) 







4 








-0.0000046598(4) 










-1 


0.0000558195 


0.000050308(7) 


0.00002159(5) 




1 





-0.0000211097 


-0.000226495(3) 


-0.0000422194 


18 


2 





0.0000249033 


0.000086519(1) 


0.0000498067 




3 








-0.0000336128(3) 







4 








-0.0000077358(2) 










-8 


0.0002681483 


-0.0206450(5) 


-0.0000340358 




1 





0.0003206506 


0.0253960(7) 


0.0006413012 


19 


2 





0.0001992266 


0.0064762(2) 


0.0003984533 




3 








-0.00199356(7) 







4 








-0.000442050(6) 










4 


-0.0004593941 


0.00110429(4) 


-0.00034842(2) 




1 





-0.0001779933 


0.00063376(6) 


-0.0003559866 


20 


2 





0.0000608900 


-0.00029572(4) 


0.0001217801 




3 








-0.000201975(4) 







4 








0.0000060939(8) 
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TABLE V: Contribution of two-loop 
diagrams to ^t]2-Ioop =1-14) 



I X 



X, 



(^) 
2A 



X. 



{I) 
3,i 




1 10 
2 



-0.0045414(8) 
0.0059942(3) 
0.0035149(4) 





1 

2 



0.0010537(2) 


-0.0012514(2) 



-15 0.0028583491 0.001177(1) 

3-7 1 -0.00492566(5) 

2 -0.001583088(1) 

2 -0.0000803010 -0.0000461(3) 

1 0.0000426116 0.00022585(9) 
8-9 2 -0.0002519813 -0.00021615(4) 

3 0.00024528(3) 

4 0.000058928(5) 




10-11 10 
2 



0.0072658787 
-0.0029032016 
-0.0010029792 




12 1 
2 





0.0003050(1) 
0.0000039204(9) 





13 1 

2 

3 



-0.00002584(1) 
0.00029039(5) 
-0.0008905(1) 
-0.000030105(2) 





14 1 

2 

3 



0.0002775(1) 
-0.0015909(4) 

0.00030041(3) 
-0.000030400(2) 



TABLE VL Contribution of two-loop 
diagrams to Sy.Lioop = 15-20) 





i . 




•^2,1 















0.0017380(8) 


15 


1 








-0.0011468(3) 




2 








-0.00090530(9) 




3 








0.000054081(9) 













-0.00016259(4) 


16 


1 








0.0007288(1) 




2 








0.0007569(1) 




3 








-0.000050771(5) 










0.0000356458 


0.00003599(7) 




1 








0.000034379(6) 


17 


2 








-0.00007614(2) 




3 








-0.000023760(4) 




4 








-0.0000109914(9) 










0.0000356457 


0.0001341(1) 




1 








-0.00022585(1) 


18 


2 








-0.00003143(2) 




3 








0.0000069846(4) 




4 








-0.0000045409(4) 













0.026645(7) 




1 








-0.04180(1) 


19 


2 








-0.00994(5) 




3 








0.0049633(8) 




4 








0.00067162(8) 













0.0010004(2) 




1 








-0.0071228(5) 


20 


2 








-0.00010774(3) 




3 








0.000068781(9) 




4 








0.0000163351(2) 
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FIG. 4: Zy^^ i^^Lf^) ~ ■^v'^^^C^lA) versus csw i^c = 3, /i = Po = 6.0). Results up to 2 

loops are shown for A'^^ = (solid line) and Nf = 2 (dashed line); one- loop results are plotted with 
a dotted line. 
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FIG. 5: Z^y'^ (QlA*) versus csw i^c = 3, /i = l/a^, /3o = 6.0). Results up to 2 loops, for the flavor 
nonsinglet operator, are shown for Nf = (solid line) and A'^^ = 2 (dashed line); 2-loop results for 
the flavor singlet operator, for Nf = 2, are plotted with a dash-dotted line; one-loop results are 
plotted with a dotted line. 
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FIG. 6: Z^y'^^ (a-^^) versus csw i^c 



3, n = Po = 6.0). Same notation as in Fig. 
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.L, Rl' 
T, 1-loop 




FIG. 7: Z^'^' {a^Jl) versus csw {N^ = 3, /u = l/a^, /3o = 6.0). Results up to 2 loops are shown for 
Nf = (solid line) and A'^^ = 2 (dashed line); one-loop results are plotted with a dotted line. 
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FIG. 



Zy'^^'^(aj^/x) versus csw {^c 



3, ^ = l/a^, /3o = 6.0). Same notation as in Fig. [71 
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FIG. 9: Z^/^^ (dotted line), Z^y^^ (solid line), z^'^^ 



Y ^u-utucu. iiiic^, ^ j^y \^ouiiu. iiiic^, ^ j\Y singlet (*-^^'^-^~^*-*^t^d line) and Zr^ 

line) up to 2 loops, versus csw {^c = 3, /i = l/oj^, A''^ = 2, /3o = 5.3). 
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(dash-dotted line) and 
line) up to 2 loops, versus csw {^c = 3, /i = l/a^, Nf = 2, /?o = 5.3). 
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